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Introduction; 


The purpose of this report ia a survey of numerical methods 
for solving integral equations. Only linear integral equations of real 
variables will be considered. There are three important types of such 
equations: 


(a) The Fredholm equations of the second k!nd 

* 4) 

$•(■*) = ^X) t h ^ K( (non-homogeneous) 

$(*) - A [ ** K < 3) 5 ( J ) Ay (homogeneous) 

(b) The Volterra equation 

Wx)* r KCx,]) Hy)<Ly 

and 

(c) The equation of the first kind 
' $(*)= 

In each of the above equations the function f(x) is to be obtained from 
the given functions g(x) and K(x,y), Under quite general conditions (a), 
and (b) possess unique solutions whereas (c) will, in general, have many 
solutions except under oertain restrictive constraints. 

Thus far, at least by comparison, say, with linear algebraic 
equations and differential equations,little has been accomplished in 
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.developing numerical methods for solving integral equations. The rupert- 
orlty ova* differential equations in expressing many physical proV ems 
has only recently hean exploited, " 

In general, it would be nicest to attack ealch equation 1; a 
method beet suited to its own structure. This puts a premium on analysis. 
The attitude taken here xrlll be to develop quite general numerical meth¬ 
ods which can most easily be treated by a computer of tho digital, type. 

On the same level of importance as the design of a solving meth¬ 
od for an equation is an appreciation of the status of errors in each 
step of the problem. Besides the usual numerical round-off and truncation 
errors, all methods of solution introduce an Important stability error 
caused by the laps® of a transcendental problem into a discrete one. In 
general, the numerical methods show excellent stability (at least where 
unique solutions to the problem exist) due chiefly to the smoothing prop- 
erties of integration. 
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g. _ The Linear Int egral Equ a tion of the F ir et .-f ind 


The equation to he considered is 

V’ 


$-M - \ KCxj<$) 


(I) 


Jo. 

with the functions f(x) and K(x,y) given and g(y) to be determined. 
Examples of this tyne of equation are: 

l 


$(*) - 


vRtt 


Z e " r ^ 


and 


'K»)» 


re ^ 




(The Fourier transform) 


(The Laplace transform) 


Often the eolution g(y) will not be unique; for, to the solution of (l) 
must be Added solutions, if they exist, of 


< 3 - <* r '^ 


> 


In both of the above mentioned examples definite inversion formulae of 
the types 

^)5(x)clx. (L, a path in the j* rlane) 

yield anslytic expressions for g(y). In general, such inversion formulae 
do not exist. The Integral equation does not always have solutions, since 
for a given K(x,y) f(x) must posses* certain properties dependent upon the 
kernel. If K(x,y) is a polynomial in x, then f^(x) must be a oolynomial in 
x. Likewise if K(a,y) = 0 for ell y and if the equation is satisfied for 
all y, then f(a) = 0. The function f(x) is not perfectly arbitrary, for 
if K(x,y) is a polynomial in x, then f(x) must be a polynomial. That is, 
we must restrict f(x) to a class of functions associated with E(x,y). 

A theorem to this effdct will be diseased later. In certain special 
cases an equation cf the first kind con be reduced to one of the second, 
kind and here f(x) can be highly arbitrary. 

As regards the uniauenes** of the solution g(y), where g(y) is 
constrained to .be continuous, the following may be *„id. If the kernel, 
K(x,y), is "definite", i.e. 

^ Kt * ) ^)vU)uiy)cU0(y >0 

for all contl nuous w(x), then the equation 

has no continuous solution w(y) except: v(y) =0. In Particular, If two 
continuous solutions g^(.y) and g ? (y) of the equation 




Engineering Hotes E=?.U -3 


Fagi; 4 


exist, then, for K(x,y) definite, & (y) ~ g (y)j i.e* the solution in 
unique. la general a given K(x,y) will not 'be definite but may hr m do 

so by the follovlng artifice. Forms 

\* KUj-t) 

Since L(x,y) is symmetric ita character!c values are reel. Fur th error e, 

(x) 1 b a characteristic function, 

\ C L (x, y) 

From which it follows that ^ O sine?® 

* [ ( L )<Mt><#') ^ 

h.no. )( rt 1^0 ° Ind.tui If the ;(x) ^ fora a complot© aet> 

^ KCx,t) 4>„(>.)Wx ±0 0:1,1, . . .) 


so that X, > O « Hence L(x,y) ie definite and the integral equation- 

u ^ x): 

Ullx): $ (i) K ( *A>ck 

has a unique continuous solution. If the function g(y) is restricted to 
be continuous f(x) will usually be required to satisfy linear conditions 


where 


of the form? 


^ UrCx') $U)<U f*o > ^ uU)jf<Cxjy)Jl}*0 


An analogy to linear algebraic equations oroves fruitful. The 
equation is analogous to the systems 

S*« t C^hKh - 


and can he obtained from ^he system in the limit as m,n oO , The prop¬ 
erties of the system depend on the relation between m and n. Hence three 
oases can occurs 
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(l) e y n implies the ’ existence of relation;* of the for.*® 


3f*J? 

X} ^X,u, “ ^ 

fJ 


for all’ y 


sjid then consistency will demand that f satisfy: 

X2lP s 


X«* =0 


In the limit these relations may taka several forms, an obvious one being;; 

b 


^ &6<) jrMdx = 0 


(2) ffi = n, A. unique set of g^a exists and is given by: 

$« = "t a * K »,^ i* 1 ' 1 ' > 

provided no relation ' 


Is satisfied., 


Kit* 

% «» K, =0 

Xt. ' (J 


(all y) 


( 3 ) « < n yields an infinite number of sets of solutions but 

particular seta may be obtained by imposing (a - m) linear conditions on 


Hence, in general,, it is too much to expect that the solution 
will be unique., An important question, particularly in numerical calcu» 
latlons, is the range in x over which f(x) io given. Decreasing the 
range any well involve making g(y) zero over u portion of the range to 
maintain uniqueness. In the particular case of the moment problem of 4t<0 
Stieltjes, f(x) is defined only over an enumerable set of values 5 x.^J 
and the side conditions Induce a unique g(y). * * ~ 


In general, then, if the conditions are satisfied which will 
make the solution unique we shall expect to obtain an expression of the 
form; 




where nay be built up of linear combinations of the type 

In two oases mentioned the inversion formula is of the highly 
desirable fora: 


‘J'.'jl- $ hl, 3, x) 


(L sone path) 
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Only rarely will this Integral «>© of the eaiao typ® as that de¬ 
fining f(x). lor a g(y) may be discontinuous at aoira y = y , and till 
yield a continuous f(x), hence the Integral for g(y) must bo inpi- cr 
that is H(y ,x) must be discontinuous or the limits bp infinite. 


In those cases where the equations of condition are not satis¬ 
fied It will then not bo possible to obtain an. exact representation of 
f(x) by an integral. Bi.it by analogy to the method of least square?, tho 

problem of solving the integral equation is essentially equivalent to that 
of minimizing an. Integral, i.e. a problem in the calculus of variations. 
The point being to obtain the best possible approximation to f(x) in the 
least squares sense by an integral of the given typo. As in the leuat 
squares proceedure we may introduce a positive weight function n(:r) which 
weights f(x) according to seme law. Then if 


or 

where 
and 
If we 


K>0 = \ b KU,^ 

i ^ 

^ |U)K(«,s1 Hrilf - j a 

4 ( 5 ) r ^ Lis.7) 

LU,<j) a y* «ca,okcm^> 

4>( S) r ^ <{U) HOgItc, 

assume that k(x,y) is "perfect", i..e. 


implies, if u(y) is continuous, that u(y) = 0. The L(s,y) is definite, 
and symmetric. Hence the solution is unique, but may not possess a con¬ 
tinuous solution g(y). This can be answered by ascertaining the existence 
or non-existence of a continuous function g(y) which minimises 


over a set of values in x. This minimum may not always exist but W may 
be made as small as desired,and thus an ^approximate representation for 
f(x) is obtained. 


A rewarding way of examining the integral equation of the first 
kind is as an operator equation; 

f(x)= dS *O.L kJ Cs> 

and then the problem is to find such that 

‘jCs'Jr 1 K f($) 
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and, In general, will not be unique. For numerical work the problem of 
uniqueness becomes largely an academic one, in any caea. The interesting 
question here ie. ’under what general conditions can a solution be found? 
how accurate will it be; and is there a quits general, easily arithmetised 
algorithm which will yield a reasonably accurate result. A general 5hoo- 
rem by Bateman proves useful, 

The equation ie: 

$60= C KL *^h L 't )< h ■ 

Either K(x,y) in symmetric or it ie symmetrized by th® operation: 

u b) - ^ - \ a K(r, x )cM £ Wxy figCy My, 

^(0 - £ b \< (s, x )KC a, y) <Lf £ )ly r £ l( y)p 

Henceforth K(x,y) will be assumed to be Byrnmetric. The theorem then is: 

(i) If K(x,y) has a conrolete set of characteristic functions, 

P (x), satisfying 

4 >„ 0 O= X„ ^ 

and 


(ii) f(x) can be expanded in an absolutely and uniformly con= 
vergent "Fourier" series 


5*(x) = Crt^OO 

Os# 

then there exists a g(y) such that 

foo- 


may be less than any arbitrary^ > 0. The proof of the theorem provides 
the algorithm for the solution. 


Proof: 

Define: 


Fb.,*.') = 



ltf-H 

Lk SCO 



HU,x> 


ftoV" 

n! 


I A 

L„ Uk) 


and 


63 1 ’5 
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The sarioa are aha. convergent for finite x» Fori 

.Lei |5Wl < 5 

IK(x,^l <- K 

They inroly 

| Lk^I ^ r $ K 

and ao each series becomes less than the exponential seriaa. 


Pick on M and write; 


W u -At 

= z£it 


Integration term by term yields oo 7*itl 

2 ^ K(x,^)F ( ^l^= a IC-0 ?_ L 

C Jz H( ^ x)dz 

- Kx)-Hfx J M> 

An expression for large M, of H(x,M) ia to be founds 

«o ' 

Uses SrC>-'> = I Crt^W 


. 7 *h 1-1 im -f a 


Therefore; 




fix 


A * l 


Form; ^ K(a^) f(^)cL^ z X 4**00 which converges. 


Also - 


l: Nx)« i ^ 4 *»w 


(I* <j l > • • • ) 
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So: 


,.. , fl -.)” Z M L" M = ft-'fz’V f Sa. fcw) 

H(x/-)= A iTT -« H Wi "6 A“ •« > 


CO 


Now 


: O I l ' I implies, for m large enough. 

A C„ 4 > r .(x) <€ 

« ' C„ 


y\i/n 


3lnce A* tl > \-* +l 


I IX* fab) I < 

n-m 'ro 


e 


l ci 


Hence: 


where 


(by formal ordering of the eigenfunctionc) 

Hfx l) - VC~0*2L ( ^2L$lnl x ) ) -f* c> 

n(.Xj2) - n 4 a? v 

I S| < 6 t~ Z/Xrrt 

- 2> «- t$. 


So rearranging: 

j'h 

HC><,Z ) - ? ^rA) £ 


v' - I 


Finally: 


so now 


1/ x 


co _ 2 /^' 

ft ( X,* ) - ^ C * <M* > £ 

i 


I f C n e' ?!/A " l 4>,Cx;| < I f c n <t>„(x)j < % 

I t f n e- fll/ ^w|<^ 


and 


for M large enough. 


Th.r.f,r. j H( X M )| < g 

Hence | ^ $(xjj < € 


and so the solution ie 


r M 


If, in addition, g(y) ie constrained to remain finite over the entire 
range of integration, more can be saido It has been s^onathat 


I ^ %h )( h ~ ■ <€ 


but not that 


lime 

6~*0 


* S (, 5> 
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If., In addition to the finiteness condition on g{y) tfe® "derived series" 

^ . -M* ( |JV» ^ . 

l~ ! Anvils convergent, $yv “ i v>\ I on (a,))) then there 

exists a g(y) defined by g(y) ~ 3 C* 6 T{y } v)&,-z RB>2: 


X 


•such that 


jU)r ^ KCx.y)^)^ 

where P(y, z) is defined by the eeriee on page [J), 


As an example cf this method of solution., consider the equation: 


e' s = ) get) aoa 5t clt 


The solution is known tc be 


The function F(t,x) is formed: 


3 ( * )s iiK^ 


oo **• 


Ft-p*>* xC*e' s ^co<$tdt - x 3 ^ts,Js ( C^s 1 .^j|e'^ J i‘'5 


2 ntf 


s e 


-t>H f^>W) x u ' + 

•fl-o H} ^ 


n ; 5 


r*>= afe-^l £ tpW 

4 \ nTc * 


cKe^*) 


X 


in ti 


Ux 


-L,-^ 

m L 

j. . -t> ( 

ifr L 


cO 


4 
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Solution la Terms of Orthogonal Functions: 

The conditions for the existence of a solution of - the equation 
of the first kind may he rather neatly stated in terfts of orthogonal 
functions. The treatment will he restricted to equations and functions 
satisfying the following conditions: 

(i) The v aid able y is real on (a,b); the variable s is real 

on the same finite interval of definition (n,b). 


(ii) The functions f(x), k(y,y) are real. 


(ill) The kernel k(x,y) is bounded and inte^rable squared 
in x and y; the function f(x) is integrable squared (L ) in (a,b). 
g(y) is L 2 in (a,b). 



Hence 


It will be convenient to say that if 




Then f(y) * g(y) 


A set of real functions 
if they satisfy: 

*b p 

<fU*> &,00d< —a, 

'a. 


f 


are ortho-normal L 2 on (a,b) 


rrm 


If the set is orthogonal and n on (a,b) it nay be made ortho- 
normal. The set is said to be complete if no function u(x) in exists 
for which y 

^ 4«(<) ^ C*) d x =• c (alII n) 

J<K, 

( f - ! V* 

If the sot*jynot oomplete, it may be completed by adjoining to it 
a eet^^)^ (13) which is called the cormlementary eet to « 

Besael 8 e inequality is: ^ 

| (Owfr« 4 s Z £ - i 

If the set f(v*GO} Is complete the equality holds. In addition, for a 
complete set, the relation 


{ f(*>gC>OoU = f 

J *- h= i U 
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is T^lidc 


implies 


I? not conmlete the equality holds if and only if 
f y 

\ %{*) <b,(*)dK*0 * ' *j 1 '- 

^ b f 6< ) cj(K) <hf~0 


\ 2 

for all g(x) on L (a,!)). This can he put another way* In order that 
f(x) he orthogonal to all g(x) which satisfy the above, it is sufficient 
that f(x) be orthogonal to the set ^ ^ complement ary to J ^ , 


The fundamental Riesz^Fiecher theorem states: 


such that 
that 


»ft.} 


cO 


? s, v < 


<PO 


then there exists and f (x) in L^(a,h) r i ’ J c‘A 


exists 


where 


C SM &0O*< 

\ t<)ck~o 


Yi-z lj 1,j • • • 

. . . , 


and if the set J <£«(*>} is complete, fCx) is unique except possibly for a 
set of measure zero. 


One point which may now be made clear is that there functions 
for which the equation (l) has no solution* For if f(x) on some (a,b), 
then for some orthogonal set on (a,b); 

C„ = £ $60 (Sbi*)<k - [ 

A„(-))r K (*,•)) & n l*)elx 

c*'i 

Art 

Is 

bo that be bounded in a necessary condition if f(x) is to be a 

a« w 

solution of (l); however it is not sufficient. But: 

x f k c f ^ k^v*. 

j k d V U 


where 

Then 

Let 


3 ng 2 Ln«»ritg Sot>3 5 ft *3.43 


?8g<5 13 


<30 t>^t Vij , < «f cohvot*v ■'■ j»t >aj f’-T' ho Trior'’' dauMHjftf; ••' •v r th" Jcr ■’ucl t . 

I* - rj *-* 

|C , • <?1 ?_ 


e> 

r- _ ’** v «— . v» 

>lttt thors ero sorie-3 2, and > ct rt which convargs, and yst for 

r-1 *3 < ^ a. ‘ 


which C, 

a} 


i-3 unbounded. Iter o:«mpl®s 


c>. 


n 


teld 


& 


S c„s r *«* -- r % 


Q.S 

i.o, both series converge. 




It io possible to reduce tho solution of th® given equation to that of 
the following problem; _ Being given a denumerable set of constants £ C J 
end a a at of functions) A (y) | „ to find a g(y) wMe h satisfies: 

c* - 1* /4*h) W5, ;*> • • • ‘ 

fhe following may bo said at once: . if the A (y) forts an ortho^corc.fil 
set on (a,b) the necessary and sufficient condition for r» solution.g.(j> 

is that 


aO 


Ic„ 


Z. • 


< <x> 


* r l 


then the solution is the limit, in tha mean, of 

N 


2 c„ *,/<»). 


a r » 


Naturally, if^A^fy)] be not closed there mist be added, to tho above, 
solution? of: 


If the 


\ c l('l)o(l-0 

* CL, 

I A n ( y>i 


are not ortho-normal, they may be .nado linearly 


independent by suppressing those which obey relations of the form. 

- 5 4 (?0 


l* I 


* V 


This creates similar equations among the c 9 at 


c* - T I 


U> 


fu~ ■ 


The remaining equations are ortho-normalized. Hence one gets’series of 
the form 


cC> 


I b, 


f (*r\ * 


h). 

3* c «) 


I 


*1 I 







( 


( 
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which must converge and the solution Is the limit ia tho mean of 

f fc, Khr 

whore {b (y)| are the derived orthc-normal set. Conversely if this nrob- 
Isa is solved, where the C are the Fourier coefficients of f (z), the 
equation is eoj.ved. 

The solutions previously given demand knowledge of the charac¬ 
teristic values and functions of the kernel or of the symmetrized defi¬ 
nite kernel constructed from it. 

It will now be show that solving the integral equation is oqul» 
valent to that of maximizing a functional equation whose maximum value, 
if attained, la the solution. The conditions that the maximum be attain¬ 
ed, then, are the conditions previously obtained which insure the solutions. 

The functional H(g) ia definod as 

where L(x,y) ia the previously defined symmetrized kernel: 

L(M): KC^*) c/t; 

fix) ia defined by: 

fa) - 

and g(x) is any L (a,b) function* The denominator of H^gjla bounded 
away from zero if L(x,y) is "closed" or if f(x) ia orthogonal to each 
member of the complimentary set of eigen functions of L(x,y)., This will 
be assumed, Define now the'Fourier*coefficients: 

C- - ^ -J60 4>U)Jk 

%- [ fa) <fco<hu 

whore tho f ( 1 )J ara tho characteristic functions of h(x,y). Then Hfr'j 

can be written as: 



«• /« v 


i h 

hx 1 \ '■ 


i 




* 


( 


( 
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Now introduce tho identity; 


u y ^ / Jjf \ l , <r 2- r \ 

l a; fk - 1 ^k + i 11 (kk - a j^ / 
> ... v «=/ <:*/ 

C-( l-» J 


hence the eeries 


A 

n=l 

<7* 

and 4? 

n-1 ^ 

«o * 





n=l 


ll < j-.f 

converges,. The substitution a. = ^C ,b f = $y / 'y. yields, the inequality 

* 11 1 


l-/ 


l M; 

Art ^ 

Now let N—9 ®° ; the right side approaches H(g) and if the series con¬ 

verges H(g) is bounded. Furthermore jfL. 'V v , converges so that 

2 

the identity can be written: <■ : 1 


JL 


00 A» 


w,) ^ J ? ( ’‘ C 4 - 3 jt i.* 

" U 2- X. A: 


Now put g, a \ f C, o Hence for the function g(x) possessing these Fourier 

1 * 1 c«_ i_ 

coefficients H(g) attains the maximum. Hence if the series ^ £. 

converges and V \ C converges, then there exists a *g(x) whose 

Kl x g 

Fourier coefficients are A C.and which maximizes the functional H{g); 
furthermore, this g(x) satisfies the integral equation (l)« 

For numerical work, one would define g^ « Aj for i f S and 

g^ * 0 for 1 > N. The N being selected as largo as necessary for a good 
approximation. The method used for the variational proceedure would be, 
perhaps the Trefftz variation of the Rayleigh Ritz method* But here again 
a knowledge of the characteristic values of L(x,y) are necessary o 


( 




H, The Method Of ?:tecpotsf:- 'Dorjcenfc? 


The aquation, to "be solved is: 

Six) t ro 

Wa fori* L (x,y)s ^ 

C jw k(x, 3 )«/x* k «* »^ x £‘ k(x^) r- - 

H'h'k): K (*,•]) K(.x,"t) 

U(<j)- ^ U t)o/t - - 

So we wish to find a g (t) such that 

f k L ( tl r(y)-0 

and the solution of this problem is the same as of (l)» 

The method of steepest descents can he ueed as an iterative method 
providing certain conditions are satisfied. It will ha assumed that L(x,y) 
is such that, for all v(y) under consideration,s 

(i) f k r(.))«lj rl-tUt V m j <ti>o 

do | [J’ L(%t) r(.t)<At^ v Jy \r x (y)ty jM>o 

(X) can he written ass 

i ft Ct/i)ir('f)d , f) >/ m ^ (’t)J't 

Now if the characteristic functions of L (y,t) form a closed set, an m can 
ho found such that (X) holds true for all v(y) which are L^(a,h). If the 
characteristic function* of L (y,t) do not form a oloeed set, a set 
is adjoined to yield the complete set+ jfty} 

Then for there to exist an m, must all v(y) he orthogonal to the complementary 




17 
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Apply the Schwarts Inequality, now, to (II) and obtain; 



- J Jf ( U^ufdydt 

tSC ^ 1 L 


j l/ 

• p 

So (II) follows if the kernel L (y,t) its L (a,b) c Ae in the case e:f alge=* 
braic equations the functional W(g) defined by 

W(^)= )‘z ) t 

<k. °~ 

is considered. Squaring 

yields the inequality? 

+ j >°. 

Schwarts 6 * Inequality appliad to the second integral and condition (T) appli¬ 
ed to the first give; 

Qn\( fty) ^y) W(y) tC tffyldy l ^(^dy >/0 

Boao„ ^ ( }>/ -Jk VV?)^ 

0 ft 

and W(g) has finite lower bound independent of g(y). 

Now, to improve and approximation, g' Vj (y), we consider ths functional; 

Yw(.,)) 

where Y, a real nurtbor, and w(y), a function, are to be detenninedo 
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Thonj |A/ ( -■ iy^ ^ ($' ( v j) + K Hf W j ( L(%t)[ J^Vt) f * (tjdt 


Collecting terms la powers of Y; ^ ^ 

W(^")= W(f7 + w“%«j 1 Ll ^ t ) X v+ ^ 


or, using an obviously convenient notations 

w(j i ' t '’)- w(^ ( '") ty+ p 1 ia/j( w w> ). 


Then., completing the squa.ro in Ys 

Wf^^lr V'/if') + ( ' /+ ^ 

- y.f w,(fX ; V) 

The largest decrease in W(g) occurs when 



?. 


Y= - */t 


(ur<"; . 


In agreement with th© result obtained in linear algebraic systems whore 
the correction vector w v is taken along the gradient of a functional W 
and turns out to be the residue, the function *r v ' (y) Is taken ass 


£ Uj,t) 


Hence, the next approximation is; 

° i u> %)= tf"h) t Yw ui (t) ^ 

- C^%) - W *"(•)) f>>)^ 

0 - /X 


4 
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Furthermore the decrease In the functional W ies 

f ^ V , 

Wlf)- ' ( 4 j_ . 

^ | b w frj (y) f *s <ri (iyt 

C f V) 

Hence the sequence/ W*g v ' is monotonic decreasing and it has a finite lower 
hound. Therefore lhe sequence converges to a unique limit. It follows that 

u,. ( W(o< w ) - ty(4 lV *°) ) =0 

y-'Xjo 1 0 u 

m I 

and the denominator eatiofios (II) so that the equation yields the result 

u>. f =o 

j * . a 

The sequence ) v' r (y) {converges to the null function. 

h«o., f 1, dt =0 

'<X 

But, hy (i) and hy the Eohvartz inequality. 

(C*( - rV'jRl Cu'.«<3 , “W-3 , ’" , '‘i)' J j] - 

>A %h, 

» m'( 

Hencet i, , , a t . , s 

|[ | [(?' | 

^ ( v \ 

Therefore the sequence, g '(y) , converges to a limit function g(y). 
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Not? if the traaisfornatlon is assumed to b® bounded i«e* eofcisifiea 



Then the " triangle 3 inequality give*? 



and, in the limit, as v—r'OO , it ia soon that the limit function g(t) 
satisfies the given equation* far 



t)|+ irfif)| £ |W%| + KUW (tj| 
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5o A Method of Br.ck Subst itution: 

The method of back substitution can he applied to the of*, 'iy Of 
integral equations of the first kind: 

a • 

$00= r LCx,y) 

’ J o. 

where L (x,y) satisfies the condition of the previous section. As. itera» 

tion algorithm, 

<j U+l Vp = f V>l 1 )+ r (L 

can be easily obtained, for 

Hence: 

where the integral operator P(y,x)0 (x) <lx is the operator(if it exists) 
that has the property 

^ b p(s,*) d *( ~ W]) 

for a p(y) in L c (a,b). In general, such an operator does not exist; and 
in every case of interest is not known; for this is equivalent to having 
solved the problem. Then replace the operator by a real scalar operator, 
' y 1 and g^ ~'(y) is definedby 

TV wj. 

This equation is analogous to that obtained in the method of 
steepest descent; however, in this case }f, being a real constant, does 
not vary from one iteration to the next. This reduces the rate of con¬ 
vergence but lessens considerably the number of calculations involved, 
Furthermore, xyith the kernel satisfying conditions I and II a ft can 
always be chosen such that the iteration method converges, i.e, 

,'X ( f (>| - f li,ho 
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for almost all. y on (a,!))- To prove this select 

- v 5(i>- 

Than, using an obvious operator formal!cm, 

<j u, (y)s V (t + <ft) fl'i )" y'i'fy)* 

Hexic© ° 

bo that 

r(t + rO vt, fh). 

Then It is necessary and sufficient for convergence of the iteration 

process that 

/4t, <- 6+ 

whei'e^olis the null operator. 

(v) 2 

Since all of the functions g (y), ▼ 3 0,1,2,... are L 

(a,b), it is sufficient to consider a function ^(y) satisfying 

<M*>= 1; 

Hence lira.(E + , v L) V 0,(y) = lira, (l + ) V 0. (y) 

V-?.o 1 *«. 1 

Since both V and f are real it is necessary and sufficient that 0 satisfy; 

O 

o< - <• 1 . 

Furthermore^ from the nature of L(x,y) it follows that, since all are 
positive, ^must be negative. Let the^be so ordered that 

— Xx * ^ 1 ^ * r 1 • 

Then, convergence is insured If O < - V< . The conditions on L(x,y) 
insure that "V, > O . 

An obvious estimate for min. = is obtained by considering 
the defining equation 

<t>, 60 = A, 


i 


\{ i t U)dx 1 Ukj'}) <t>,(y)cii , 

Jo- . < 

I 

},(ft ! cCk; i^jt x ) Vv £* 6, (*)If 

O- «t 

or, using an accepted notations 

1 > -L 

/ || LU, ■j)#. 

Hence e possible ^ is ^ 

Furthermore, even if L(x,y) is not "closed" but is "definite 1 , the above 
estimate will insure convergence. 
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Then; 


^ <k (^’<4 


By Schwarta * inequality; 
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Solution by Quadrature? 



A rather obvious approach is t.o reduce- the integral c.v on 
to a set of linear algebraic equations, whose solution will presumably 
give an approximation to the function g(y). This may he done in e*y of 

several ways. 


\*£ 


Firstly, the Integral may be approximated over some sot of values 
by a suitable quadrature formulated. As an example, use T.ey be 


made of the Gauge quadrature formula. The range (a,h) is transformed to 
(-1, l) and the integral is approximated by; 

£ KIm> 3‘1>^ - £ CK(x /3 .)g(y A ) 


I 


where the values o(^ ^ are real constants which depend on the number of 
points chosen. The set of values^ =n are the n solutions cf P^(y) = 0, 

feh 

where P (y) is the n ‘ Legendre polynomial. These seroea ara real, fiie<= 
tinct, and lie in the range (-1,1). They have been tabulated for fairly 
large values of n. By thie double choice of weight factors and ordinates 
the Integral, evaluated by n points, is exact if the functions integrated 
are polynomials of degree 2n - 1 or loss. For small n this is a decided 
improvement over the Lagrangian inteipolation polynomial. The given func¬ 
tion f(x) is then aval tinted at the Get of points £ s^ 31 * with each x^=y 

il 23 ! 

i=l,2,...n and the system of equations? 

$C*p= fat; 1 "' K(x Jy y.) 5(3.) ^ 

ill 

is solved for g(y^) (i = l,2,...,n) the answer being obtained a« a set 

of value? ^ g(y t ) ^ « The integral equation may be used to afford 

some idea of the error. Once tho approximations are Introduced, there 
is little point in debating over uniqueness of solution. Obviously this 
method may well give a solution to the algebraic porblem where there is 
none to the integral equation. However, if no relation of the form 

.A 

c - * 

ie satisfied, this method should give a reasonable approximation com¬ 
mensurate with the number of points used. In actuality if the rangec 
In x and y are different, a set of linear equations with a singular 
matrix will be obtained providing the intervals in x and y axe 
different. 
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In additicxi to the method mentioned any of the iterative so ham as 
for- calving simultaneous equations may he used to advantage. 

However the important faster of solution stability mot b- kept 
i n mind. The reduction to a system of algebraic equations involves one 
approximate quadrature and the remaining calculations are algebraic. An 
iterative scheme like the steepest descent or back-substitution methods 
Involves repeated integrations each of which Is approximate but essential” 
ly no solution of a set of linear aquations. 

It is a matter of extreme complexity to make a precise "error*' 
analysis when an integral equation is solved by numerical methods; however 
the need for such an analysis is pressing. Of course once a solution Is 
obtained the error is immediately calculable from the equation itself. 

An a priori estimation of error, i.e. of the needed fineness of the in= 
tegrntion scheme, or of the degree of approximating polynomiels, however, 
is usually a problem of the R;>ae order of magnitude as the solution of 
the given equation. Since iteration methods tend to forestall the un¬ 
limited accrueing of round-off and truncation errors which occur in each 
calculation, It appears that an Iterative method would show more stability 

The functions can be maintained as tables of functional values 
for given ordinates or as tables of coefficients of interpolation poly¬ 
nomials, in which case presumably the quadratures can be performed exactly 
except for round-off error. Certainly, the ca.se of the kernel, a function 
two variables, if 100 points were required, HP values need be stored, 
whereas an interpolation polynomial of PP degree would require as much 
storage for its co-efficients. 


Rnt ii^3rirg l>c%m JUlUj 


f 


Vfa vieh to ol>t?*ln & numerical solution to «&n integral iiqn^tion 
of th& tvne 

f>cO 


$u)* ) . kcx,«i)^w 3 

** CO 


This method will bo limited to those earticailnr cates where k(r,y) is 
of the tyr-e k(x*y)» k(x+y) or k(xy). 2y a simple variable transformation 
each of these can b® transformed Into k(x-y)« So we wish to solve 

r <70 

$W= j 

Formally we can obtain a solution in the following way . 

Formally; ^ 

F(u) = V 9-l*)& f ' UK cfx 


'JTtt 




1 'i£r re u V* C 

vrjvr $*° j KU-y)e''**M. 


S'-&& 


let x * t + y and obtain 

r , . p’® , . . A''’ 0 4h( tf») 

Fix)' ir=r i \ K(*>e c di 

' /V ^d© -»oO 


Therefore; 


^ X» ^ e ‘ v " L ^ t;e ^ 


f{w) "ST (Jtff (|N> «(ut) 


where G(u) and E(u) are the Fourier transforms respectively, of g(y) and 
k(t). 

Finally obtain: 

/*, x J- F(w) 

w* u - wr KU) 

and our desired result is 




j. rpfx) -" 1*4 

itf J l/l,A ' 


^ KM 


w 


The recvl»t io purely formal. The convergence of th® ismropor integrals 


and the validity of interchange 

of in 

tegratioa ordsrs must ba in'. 

The fcllowing theorem, 
Tv is 2 

however, 

win 

')« of use in many numerical 

Theorem: If: (i) 

iU) 

£ 

L l 

(ii) 

KC<) 

e 

t 1 ( 

Then if and only if 




(ill) 

F(v<.) 

KU) 

£ 

l l ( - ob y oC ) 

is ^I'i) £ (. l ( 

*eP, oo) 




Ve'nre specifically interacted in a numerical solution when 
the given function f(x) hae certain properties: 


(i) It is defined, either through lack of complete information 
or for convenience, at hut a finite number of points* 

(ii) These points, except perhaps for an assumed constant value 
at oo , are distributed over a finite portion of the real axis* 

(ill) It will usually be the case that the vsiue3 . tabulated 
over the set cf points^ x A \ will not be ffj^) but f(x^) + g (:c t ) or even 

more likely f(x.) + g (x), where g (x) Is an error function* That is, 

the simplest case is where the error at e,ny noint is functionally 
defined by 

eu,)= 


Usi\ally, through having interpolated by one means of another to obtain 
each f(x.), the error at any point x , will depend on all points, i,e.. 

Naturally wo are going to assume that a unique solution to the 
equation exists* And, clue to the incompleteness with which f(x) is 
given, we must eroect to get at best an approximation to g(y). 


The Numerical Solution 


fw= u# T f aic iw <(» 

liiirFM- \ * r 


(ft) 




( 



C ; 
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Pass' 2 ••? 


The cut-off at x 


la the integration iatar^al and the 


subsequent dropping of l(x ya) in equivalent to folding the tr&v ■ form 
P(*r) egoInst tho transform of the unit atop function: 


m(iO= ) ? 

Lb 


X > yL 


r *> , . x < - 

i.e„, ( F( t) sin ^ °^Sgfe} (u-t )dt. 
^ - 


; •' * ro 
‘X© i X - *0 
X <• -X 0 

__ The praciae affect this ha3 on 

— co Xe 

tho solution ie difficult to estimate. Hones wo have here a first li¬ 
mitation. It is necessary that a finite r exist such that I(y v.) 

O 6 

will bo sufficiently small. Since F(u) exists this would appear trivial> 
but if f(x) is a tabulated functicn. or an experimentally measured quan¬ 
tity, knowledge out to x q , at least, is required. It is possible then 

then that the range of definition of f(x) may be so limited that l(x o }U) 

will be of such a size as to a render a good solution impossible. This 
difficulty will be discussed later. For the moment we will assume I(x ,u) 
to be sufficiently small. 

The evaluation o? the transforms requires that there be avail¬ 
able tables of sines and cosines or that they may be conveniently manu¬ 
factured in the machine during the procop^ of solution. We assume sets 
of trigonometric functions will bo available. 

Further, on the x scale, we assume, since this will usually be 
moet convenient, the values of f(x) to be spaced equidistantly- then lot 


and let V k 


ITT A.* hu/* 


. • • h. 




iffir 


0 ; °i k j x > • ■ • 

Then K — X . We have: 


pu > 3 5,7 h-xOw^) 




We can alws/i- make dh = 1. Using, in affect the Eulei*-MacLaurin summation 

formula; , 

v hm j* • I 

Ao i 1 1— _ _ itr c idle 




e, 5 : w„,< 

Similarly, for the kernel, 

Jii? KU) 


a A 
b 




X* ^TV, , -7 

* c 7 e 


-l 


may. 
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J'orm the quotient: 

(t 


fcU.K 

z c « 


urihw 




_ —rr - .8; : j = 0, 0h 

Ka4K _ -(JTt.W 

l l w*>;u - 

jt )C I { i , 

which ia to represent ? r Ui ) at values of w corresponding to u 

That meana that for 0 &. w £ v « u runs over C4s£ ___ dy ' h\t « 

max' v ,¥ MdK 

Heat wo form 




Aur 


<3W)= ~:\C,Mc'^cU= b C Zi*> e'"f*cfK 

( 1&J 1 ITTT -L k7*J 


co 


oo 


— cfc 


irr 






riwi e -‘f M . 


UiV XU) 


t** T (<w$). 


-u KM 


The Euler-MaeLaurin Integration scheme gives: 

. . v u,*** o r p( U ) 1 -itCmz 

” s'’V.* »=. 5 ,., [ a£L.< 


end again let A m = 1. The variable rangaB are as foll<n»*8: 


(i) 


x: 


0 to x. 


(ii) h: 0 to H(» h ); A h * 1 

max 

(iii) v: 0 to Vf(a v ) (apy 0*250 of a cycle) 

Jftax 

hence from H vpluoa of the integrand, we get 


W 

values of ?(u)o 

(iv) u: 0 to 


ZTT HW" 


X 


W 

Aw‘ 

subset ia evaluating g(y) 


3 


of the H values of [ EJii^ j we.may use all, or an equidistant 

L. KCw) -*u—>«/ 
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Hone® while (iv) holds, 

(v) ra: 0 to M in H jumps: A » = 1< 

w 

&od ----- ss Ml; 1 an integer. 

(vi) 2 : 0 io Z( = 2 ) 

max 

r p/u)l 

hence from M values off U-- we obtain 

L klu) -1. 


Ktw) o u-^yT 


Z/ A * values of g(y) 


y-^2 


for which: 


(vii) y 


irr/^l-Ly - /1Z / 

: 0 to ,rrrnr Tjy/ 0 


z«HW 


Having obtained g(y) , the problem is essentially solved. 

More may ba said about the actual evaluation of the transforms. 
We can always choose our zero on the x axis so aa to maintain the follow¬ 
ing symmetry: 

F(*)- ( j=j. ^ j [fait $(-*)] Ort** + A.[f(x)- ? c/k 

and likewise for K(u). Hoth F(u) and K(u) are evaluated over the same 
sets of -ooints. And 

, j_ » 1 . x ( P~« tcF,X *(.-<■ Kt ) 

1 ^ G Kc'U) +. 


Since we Integrate over (- (fi , oo ) all odd terms vanish,. Hence we must 
furnish: 


(1) 

*c ( “> ♦ 

(ii) 

F K 

c c 

(ill) 

--IT K 

IB S 


C 


and 


g(y) ovor the range 0 to y v will be the (cosine ■+• cine) trannform of 
G-Cu), while being th© (eocene - sine) trrnsforo of 0(u) over ~y to 0. 


mao: 







( , 
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Cut-off error s; I n genera 1 ., 


It is very difficult *0 osaee 
ih error In the solution of this problem. The way they arise : at 
once evident 'though their accumulative effect on the solution 1 


».U C 

io.iilt to determine. 'Hie errors are tho usual ones.’ in a numeric 1 cal¬ 
culation. Wo may Hat thorn hero. They are: 


(i) Tho function f(s), either through oxpadloncy, or in 
fact, is incompletely defined an, as mentioned before, .will never, ox- 
copt possibly fer a constant value, b© given over the entire rangs(*» (P, CP) 

(ii) In addition, each value, f(x.) will be contaminated by 

1 

errors.. These will be, at least, truncation errors. 

(iii) The infinite integrals met always, unless they can 
bo evaluated exactly by pane analytic procedure, be cut off at some finite 
value. This is equivalent, as mentioned, before, to tho folding of the 
precise transform against the stop function t id UK . if x is the 

’ u.y 0 

cut-off value. This means that the frequency spectrum of f(::) ia warped 
by thie factor. What is more important, though, the warping cannot be 
removed form g(y) by any unfolding. 

However, we car. make gore remmptlona about f (x) for large X 
and develop an error formula for the cut-off. We observe that we have 
two cut-of" integrals to evaluate: 

u.) 


<b) 


S fMe^dx 

J X. 


As x —9 oo , those errors —? 0; and the high frequency components of 
F(u) should be largely transformed at any reasonably large x and similar* 
ly for the cut-off error in the evaluation of K(u) 

The integration scheme 

Tho evaluation of F(u) and <l(u) is to be done numerically and 
in thee© cases where it cannot easily be evaluated analytically K(u) will 
be also so treated, 


Since ve are evaluating special tyoea of integrals, i.a, Fourier 
transforms, it is best that we us© an Integration scheme which allows us 
to make most use of the duplicative properties of sines and cosines. 






( 


t>3 h S 

J'nginsuring Ho tc s 


rage 


To make "best u<se of the computer's limited storage capacity tfe must use 
a process that forms the Integration components from p. minimum storage 
list of sines and cosines. Shr the moment vo vill stand pat on the Buler- 
MrcLaurin integration scheme-.Asing eq.uidista.ntly spaced points. It is; 

7j Trt-» 7) 

\ j +• f(n) - j f (r) S, ('t)dt 

* C~ m "i ”'i 


n-2. 


vrher® 


<*> 


-g I, 




j (/.A 1 TTJ1 t 

^r>rr 


V . r 9 are to evaluate an integral of the form; 

i °5-60e tux d* - ^ ^ 


x (T t h w 


dfi 




H * *. * h ^y 


The sum formula yields; 

*V| H 




] 


n-- ( J Q 


An analysis of the error tern is desirable for it develops certain facts 
about the relation of the h end v intervale. Naturally as the number of 
intervals —oo , the function S(h) —t* 0. The functional relation¬ 

ship between the error and the interval length resides, of course in 
S(h)» We will assume that |f(h)l < M, M a finite constant, over(p,H). 
This is not unreasonable. We consider the error for W =* o» 

rtp ( f ^S<Min-nh)dh 

0 i 

: M f "(n- - CO )dh~ 0 

J o 

([h^ie the least integer in (0,1), (letc) 

For W a fth, til = integer 

A / j\ M f H urtrwh . , h . MH 

£W- Muir)ui«\ e ^ tT ™ hdh ' wrr 

n=, 




( 


( 


*ing 


O.gO 




For th© error to to zero. except* at \f «00 , it is cecostary that a 

over (0,H). Thin it cannot be. 


•>.R >■ kwf 


be orthogonal to 

a* 


y s^tj 
~v5f 


Hence it would appear that, sines ’rf « 0 is a va'ue wo-use, we should 
though tho peak about W » l is narrow,'Unit W to lie certainly in the 
range 0 *» W £ C < 1, to minimize the error. This of course in a con¬ 
sequence of the fact that we are fortunate to be dealing with periodic 
functions, q'H.j hw 9 super-iinposed uron, as we have assumed, a function c 
containing practically no frequencies above some finite cut-off point.. 
This concludes the analytic discussion of the errors exclusive of these 
occurring from round-off, truncation, or in tho function f(x) itself* 
These latter can probably be treated statistically and this may actually 
be done later* 


Convergence Factors 


Before going on to the more numerical treatment we must men¬ 
tion some important considerations concerning convergence. We assume 
that F(u) and K(u) exist and are finite. Hence F(u) must have a higher 
order aero at infinity than K(u) to insiire convergence of 



But it is precisely those features which we lack in any approximate 
solution which make this extremely difficult to achieve Every error 
that eccure in the procedure has the effect of adding into F(u), 

Fourier components which extend the range in u on which F(u) maintains 
appreciable value. It is only through a sacrifice in "resolution 11 of 
g(y) that wo will be able to obtain a convergent numerical solution. 
Nevertheless, admitting tho necessity, and even Justice, of such a re- 
solution depletion wo can obtain solutions which are quite good. The 
extent of this "smudging out" of g(y) will depend, then, on the nature 
of g(y), the extent of definition of f(x), and, of course, the accuracy 
we maintain in the numerical computations. A datailed example will be 
shown later where the function g(y) is a sum of J - functions. Even 
though its transform does not converge we are able to obtain a partial.- 
ly resolved solution; indeed in this vorst possible case we do obtain a 
resolution as fine as desirable by merely extracting an increasing num¬ 
ber of values of f(x) and refining our integration scheme, accordingly. 
Hence to obtain a g(y), we must usually expect to be forced to introduce 
a convergence factor. A convenient one, since up to cona^ent amplitude 
and half-width factors it transforms into itself, is e" 3 ^ . Another is 

S lav ^ a * and still a third is X 1 . As in the formation of 

ICAX)*' 


i 
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. 3 

'Bn&zeortz 

tho first 

transforms f( 


PfiCV-V 


If < 1- s 


■factor givas as a solution at the 6"; f. r :*.c of resolution, For, 
cana t oalcnl at © 


vre 


C (v.'i - ~ 

m, A i '**• m 


i r 


o c * 


/-Tl 


Fiu) /C ->^ , 
77 T v c ' cite 
7* Kivt) 


We introduce the factor T(u) and get 

¥ 


a h) 

1 J 


x. (' 4< V 1 ^ > -r(*\ ^ c t u ^/ u 
LIT 3 77r: x ; ° 


■ tir 

. f 


-'cO Kt*) 




- \ § (. ■Jj - $) -its) 4 5 > 

. A o® 

In short, our solution is M folded against” the transform of T(u)= The 
step function yields a g* (y) which has oscillations of decreasing amp¬ 
litudes and obscures much of the result. The function e“a.u transforms 


into 


,> - % 


, up to an amplitude factor. This means that tho 

' FI*) 

greater the value ia u at which well behaved, the narrower 

will be the spread by which each point will be weighted in tho £*( 7 ) 

This Is almost poetic in its justice? . For In just this way does the 
incomplete definition of f(x) affect the result. 


A s pecific example 

Consider tho Laplace integral!, equation 



Kf )= ) 6 



The transformation 



- d _ X 

V-< > f I£ 


yi eldfj 

^ •* A) x - «d 

e < $(»*)= \ e 


or 

- c*> 

i* iw* S. « e 

f • 

We forms 




¥> 


kw- i- r e '*V«‘^t 

\iW " 


V 
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! :!t> ?ir.e 




0.^3 


and ’chair quotient to gives 


aj , . / K 5 r/ u i o M U 

n (uj\« — \ f ' cIh, 

i «' W -V l««) 


With thie kernel, K(u) can he evaluated analytically find is easily soon 
tc ho the Euler integral forJ- f*(l «*> iu) . Likewise a necessary and 

<K<\ ' ° 

sufficient condition for the behavior of f(x) for large x(i.e. x ) can 
be stated. With respect to the Wo cut off integrals it is no follows2 
The cut off integrals are: 

re'*?(-x)e' l “^ 


(a) 

(V) 


o 

.*•* 


L eMu) 


>. L U /- 


ck. 


Then for convergence,:, it is both necessary and sufficient that, for (a), 
for large x, f(~x), bo Oie * 1 *), n < 1 or f(x) be 0 (e“ M ), n< 1 <; and 
that for (b), for la.rge x t f(x) be oCe 1 ^)^ n < - 1 . For both these cases 
an error formula can be simply obtained by a single integration by parti 
They ares (a-O X 

(a) j iV'Olf" , )e-“V<|‘| 1 -i«(f ’>| 


j ^ e-*6(e‘*)e-'“% 


and (b) yields similarly; 


If e * f(K) j £ Qj 

*♦ jy\ 


fiiiLK, 6( 

\n-J! J *o 

Q( e (v '' ,)Xc ) 

I » * t j "" I t* I • 

C nt OXe 


(V»>X 4 


- Ut X 


)f ax 


|>i +! { - |u( 

A more specific result holds for the cut-off in forming K(u), 
W® observe that these Wo formulas yield an exact expression for the 
error in the two cut-off’s. They are; 

$ e'^e* e lU cH - -fly 

t 


and 


h X, 


oO / > 1 K 


• where r,M referee tc tho in- 

complete gamma function of x with Integration limit j. 
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we observe that? 


. f* ( \J. 1 

•Having ferried the quotient -v.,..- ■*••'/■ 

\'(t+vn) 

j rt i+i.u^| = \S. |s ^x TU _j 

so that K(u) diminishes quite rapidly. Assuming that a coavo.-. »ai;a© 
factor is necessary let us introduce e"* aud as caa® Our choice of "a” 

f*» f ^ ^ 

Is governed by the behaviour of "7 — v aa u increases. The necessary 

s<WW 

existence of m\dfflmpod Fourier terms in ?(u) implies that for coma 
value of u, the quotient will begin to grow without bound. *•& select 
an n r" such that ( 9 ~A.tc v will be essentially zero for all 

v<wr 

u beyond this critical value. Hones v/e have? 

" f?u > 


•" K, s 

H 1, P = 


1 - v 

a J 


-oTi 


_.) 

r(iuO 


e e J 


cJu 


1 ffv f ?‘rt> 


and we eea at once tho way in which a snail M a* (i.e. a largo critical 
"tt") yields a good g*(y). Here then is the crux of the resolution 
problem. It Is all bound ut> in the constant "a n . The* size of "a ' 1 will 
then depend u-ian just how well we can damp out the M nol sj 1 in the evalu®» 
ation of F(u), That it can be made to be as near zero ao we wish is a 
consequence of tho uniqueness of the solution. 


J 


i Inal; 


E»K3 


s® jr 


' > 1*b ! 3 t j Set 

’ 


Shs equation to be const 

p (9 

ur( i ,0 - 3 (./•) + a \ 

U flL -' 


doled Is* 

KU,i)urlt) 



She functions v(x), f(x), aad h(x,t) are assumed to he real 

functions of the real variables x s y. As In the prscoeding section,, if 
&(x,t) is not symmetrical it may be made to yield a symmetric kernel* 
The above equation may be written in symbolic forms 

wCx)z ^(x) -f }T w(>0. 


01 ](.*)-- (e- >tw*> (*) 

where £ i» the identity operator and T is the integral operator* Ob¬ 
viously the equation is solved if (S - AT)"^ can be obtained* In general, 
it will be shown that such an inverse exists, subject, to certain condi¬ 
tions, as an infinite series in powers of T(l*e s as sums of iterated in¬ 
tegrals) , 


Left multiplication of both sides of (*) by (S -- AT)’ 4 (the trans¬ 
pose) yeilds 

tU*)= ( E* M) j(*) (t) 

where L is a syranetrie Integral operator and hence has only real eigen¬ 
values. The equation to be considered then ie either (i*-) or (t) * For 
theeo eqtiations the following existence theorem is fundamental 

If the known functions g(x) and K(x,t) are integr^-ble sqtiared 
over (a y b), then is w(x) integrable sciuared over (a,b). Furthermore 
if A is such that the homogeneous equation 

4 >U) S aT$(*) 

has only the solution 0{x) = 0, then there exists a unique solution 
w(x) * If, however, the homogeneous equatio. lias ns independant solutions 
then there exists a w(x) if and only if g(x) is orthogonal to each of the 
a solutions of 

A //( y )- AT'fCx) 

The solution w(x) is not unique for there can be adjoined to it any linear 
combination of the m solutions of the homogeneous equation* 

2 * The Simple Iteration Procedure 

Va form, from: 

uU) = H* ) 4 A [ u 

• la. 

the iteration procedure; 


U) 


u‘"Vx)* f 60 + xTu c "‘' i 




then 


v'&jre 


u, 6 o - ft*) h a T k*)+ a r~Rx; t . . . 

A *~« ^ 


K«* fr* u b u) 


let J v (sijjtfcU,, then 

|R.UU’llui» r * 


a necessary and sufficient condition for convergence is that lime II R (I s 0 


00 that: 


V, ^ 00 


im I ri«T 1 -*»o ; 

n-^ou 

ft sufficient condition is £ of course* 


I xi im < 


ors 


I MULII < I 

We define the nsoria of T by the following equations: 

II Lll - ( ^ I 


a. <L 


and 


^ LCx, v^) a T 1 


We define; 

I ITi*0<) lt ^ lf jn 

II uU)|l ” l®** the least number for which the 

inequality is true. It follows that 

IITull 2 £ IILII 2 IUM| ! 

from the Schwart* ineouality.. Hence;’ 

IITil * lUH 

Furthermore 

HT k u(oU IITI l|T K "uWll, 

Hence j if ^ 

m- T T’ 


r^i 


we have 
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The solution, so obtained is unique. A major problem; then, is how to 
make use of an i-tei-ativs prooesB when, say, |! X HI j 

•bet iiXL(x,y) |j be bounded. Then an algebraic procedure can 
bs adjoined to the iteration procedure to provide fe solution when 
|jAL(x t y)ij>| o It is always possible to find a set of functions A,(j;:) s 
such that; 

I! L(*,y)- X AWg^llO 

or *" ‘ 

L(x,'j) = t'AO'f'Bji) +Mix,</) 

A- | 

where 51 is chosen so that || M(x,y)j|<.le Hence M(x,y) has an inverse, J, 
consisting of convergent infinite series of integral operators Then* 

U.(y)rTifx) + T(7A t U) ) C 0*tf) 

1*1 0 


or, using an obvious notation; ^ 

( 8 J (sO > uC<)') - (8^00 , T f60) -b X (•?/*>, 74 1 GO ) ( 


(PjU>,u(*>)« £ uU)cAk 


j - 




The solution this set of linear equations for (B.(x), tf(x) ) when sub¬ 
stituted in the equation for u(x) yields the desired solution u(x) , 

This method, in essence, can be regarded ns a double Iteration scheme, 
if the algebraic equations are so solved, 
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Tbs Method o f Steepest, ^o soant 
Tie equation 

$60 = ^tv> + .x ^ k(*,y)ft'f')dy • 

in operator form is 

(£-XT)^(*> =^(x) 

and T will be assumed to be a positive operator* Then, it - follows that 
the inequalities 

o< yv» cHpC^V)^ 

are satisfied, the first from the positive nature of the operator, the 
second from the bounded nature of the operator and Schwarts's Inequality n 

The functional v(f) is defined by 

W(f)~ { /z 5 K*) V - A fM ^ ) $r(y)ely + j a f <*> f)C* ) ^ * 

Squaring the expression 

-X $C*>^KCx,*?) f (*)%<* Wx. 

yields the inequality 

-A Kdx, j ^ Jk)g(*>cUy 

from vhich follows 

5m( S** SVx)dx) W(^) 4- ^ *0. 

Therefore w(f) has a finite lower bound, independent of f(x), given by 

w(/)» - L" ^ t>) o>x 

(v) * r>1 ^vS-l) 

If f (x) is an approximation, then an improvement, f ' ' (x) 

is obtained by considering 

W(^‘ ] ) - W( ’f l "c«>+ 

(v) 

where a real number, and r s '(x) are to be determined. Then 

W (W If (v) ) + ^2 C 1 A " >(x) 

+ *[\u) /f'fadx +• <**T $n w, <x) - >n u; ^)£ 

Jce J °“ 
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or, in shorter notation: 

W(f" ,: )= W-f'V f <*©, -r Of*®* . 

Completing tha square in ot, s 

w(-T’J * © u (<* + J4^f- 14(f /. 

w(f) decreases most rapidly when 

«. - - 'k 


as before, the function. 


& 

, , & 

r (v !x) 


is taken as the residue, i„e< 




The next approximation is 


$ l ’ ,n (x) * P"U) r«n M (x) 


and this is the iteration algorithm,- The sequence, > w(f being 

monotonic decreasing, has a finite lower hound and so converges to a 
unique limit, from which it follows that 

im [xr'*")-w(r')]-° 

y*>va 

which hy a previous inequality implies 


-U- CVW'J*** 

1 / —wa j- 


V —>c« 


from which it follows that £ f^‘ (x)} converges to the solution f(x)j 
almo t eve rywhe re „ 


( 
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J Pj>- ‘ J C- 


10 e l ‘ h a Method of , Sack^Sn b*-tl b utto n 
The equation to be conaidersd. is; 

Hr)s^(^) t A j^KCx,^) ffvMy 

lito operator in question >nay he made ^definite “by reduction to« 

-J(^) ~ 1 + * p KC»C/1 ^^Mk J 


f 




}(^)r Ulv,) 4 #*)«** 

The iteration scheme 1st 

(✓>' V . . . {*! 


and 


whore 


f (l (<})+•* ft 

nTb)* 4 ""(<*) - A C LU '^ r'V'Mx -;/<■*) 

and q[ is to be chosen so that linu | r^ V ^(y){ = 0 

V** 

Than* , 

A <w ", ¥ >. f "" V- * I **>•>* - 7^ 

; f'"iy + (AA"'(>)) - A [Afx.^D^O 

k 

n. utu (y)i A""h) + « Ao^lU.-P n '■''C x )c/?c 

or, in operator notations 

a <' , -" , (v) * ' (f + *ce -ATJ a‘"("))- 

Since the operators are bounded and commits this becomes s 

d'^'^V [ <?+ *(p-at).3‘Vy 9 j 

• ( 1*1 \ 

find o( must now be chosen so that Ilia r ' ‘ (y) — 0, Since the functions 

2 K**> 

considered are L (a,,b) and T is a bounded "positive 15 operator* it will be 
assumed that 

0< /y\ C ^°(> L 
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•?£.ec 


■‘i3 


P 2, 


tor an y f(x) la I/'(a,b)» Prcm the above it follows; that 

|! n""'}( 1 )|| £ || 9 ■* Ue->~T)\\ V 
(y) il = o 




Senes lira 0 fj 

CCi 

if ancfonly if 

(| t + < l 


or 


^ f < o'vX iU,^)U*f)<fy < L 


Then, 

| I +• (A “ oC Xrti j < I 

yield* a relation for <X „ At this point it might be mentioned that the 
convergence criterion is the same as for Newton" * method.. witho( » - y 
Hence, BUbject to thie change in sign >: the two methods converge or di¬ 
verge together,. Thie method also fails, then* for^ = X* > a charac- 
teri8tic value of the homogeneous equations Newton 1 ’s method will, in 
general, however, converge more rapidly. 
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lie The Solution by Kjsrton* a Method 


The equation 

^(.y, )r <^L'< ) + X ^ •' :rt"t )c (“t~ 

can be represented in operator forties 

( c- XT) gC*0 

defined on the interval (a,b), where E Is the identity operator and T 
is a poeitlvo integral operator* An operator (E ~ A T)" ? exists and 
in the case of it AK(x»t) !! < 1 is rapreentabla as an infinite series in 
iterated integrals operating on g(x). Newton's method n ay be employed 
to obtain (X * AT)“~ » in those cases where 8AX(x,t)ft £1. 


By analogy to the Newton algorithm for obtaining reciprocals 
of numbers, the equation 


pw-nK _ c e - > t3 


i a formed where the operator 



ia ouch that 


li. p<" = (f-at)-'. 

k(o) y- ' 

Considerations of commutativity may |>Ov ignored if p - o , a scalar. 
Then p , and hence, by induction, p' '(v =» 1,2,,..) commute with the 
integral operator XT. This is quite obviously a oufficient, and oven 
convenient, but not necessary constraint on p' ' (v = 0, 

By applying tha recurrence relation it is seen that 

e- (£• (c- - ( 6 -*T)p (, ’) ! '' 


Hence 

lira f ?t/l ~ C £ ' A T) 

if and only if 

lim II 9 - (<r‘ XT ) II Z r 0 

V -*♦<*> '' 

A sufficient condition that this be eo is that 

|| £- (f-vrip"" IK 1 
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■ y 

T.y a proper choice of p' ' thi* cun always he arranged. Hence Now ton ? s 
.:-atho4 »dl always yield a solution by • "urely iterative scheme ' v f •••rid¬ 
ing of course the operator T ie bounded, which will be assumed s.‘. all 
times). Obviously if A = A j ie a characteristic v£Lue of the equation, 

the inequality above ia not satisfied. Likewise such a A violates a 
necessary condiMpa f i.e. that all characteristic values of the operator 
E «> (IS « \ T) p'° ; lie outside the unit circle in the complex plans for 
with a A = A j the operator has a characteristic value on the unit 
circle. 
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13. Soluti on by A pprogj.uat a Quadrature 


The Fredholm eeuation 


$G0 s j(*) + *£ KC>',y 

nay be regarded (aad was so treated originally by Fradhola) as the 
limiting ease of n simultaneous linear algabi’aic equations 

K*;)-- ^JCj+(Vrf-r, *. 3, yV 


as a -^oO„ If >is not a characteristic value of the homogeneous equa¬ 
tion, the above, for a large enough, may be assumed to yaeld a solution 
f(x^) i= 1,2, cc, a, which closely approximates f(x) at least at n points. 
Since a. soluticn can only be obtained for a finite number of points, la 
any case, when numerical methods are employed the above appears to be an 
attractive way to dispose of the Fredholm equation by reducing the trans¬ 
cendental to an algebraic problem c 


The above equation represents a solution at "n M distinct points 
equidlstantly spaced* However, it involves the most primitive quadrature 
method* A linear quadrature scheme which best approximates the integral 
is desired. One method is to approximate the functions involved by in¬ 
terpolation polynomials j, and prevail upon an exact integration to give 
the algebraic system. Lagrsngian or orthogonal polynomials are most 
commonly employed. Another method involves the use of non equidlstantly 
spaced points and the concept of Gaussian mechanical quadrature. Firstly 
the interval (a,b) is mapped into (-1,1). The integral is approximated 

5 r X c K(x,^p 

(n) 

with a given fixed n. ThetL' ' are weights which vary with J and each 
n. The points ,y( n ) . l,2,^..n are the n real, distinct, zeroes of the 

nth Legendre polynomial and lie in (-1,1) • Tables of the y\ n ' are avail¬ 
able for n$ 10 and, in time, will be no doubt available for larger n, 

Such a quadrature yields an exact evaluation of the integral if the in- 
tengrand is a polynomial of degree 2a - 1 by the use of only n points 0 


The method may be generalized to treat any interval providing 
the points and weights are chosen to depend upon that polynomial ortho¬ 
gonal over (a,b) with respect to a unit weight function. The more^oi^ 
thogonal functions may bo used providing the proper weight functions are 
introduced. For example, on (-1,1) the Tehebichef polynomials may be 
used if the integral is written ass „ _ _____ 

c' t i r Y, 

-i ij&p' t=« 

Here the j3^ are identical for all i and a given n, though the yf°^ are 
rot equidistantly spaced being the n zeroes of the Tehebichef polynomial 
T (y) . Likewise use may be made of the Laguerre polynomials over (o,os) 
w?th weight function e" 3C ^Vnd Hermite polynomials over (-co, *o ) with 
vei ght e" x ^1 % . 
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